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Nomenclature 


A,A',A",c1,c2,c1t, 

c2T’cT’cel’ce2’°e 


u.v.w 


constants  In  the  model 


drag  coefficient 

body  diameter 

Froude  number  (=U  /D  N) 
o 

acceleration  due  to  gravity 
acceleration  due  to  gravity  vector 
wake  height 

turbulent  energy  per  unit  volume 
Brunt  - Vaisala  frequency 
mean  pressure 

pressure  deviation  from  hydrostatic  state 

Richardson  number  (R-1  = r /$  ag/V 

o o 

Richardson  number  (=  HV/k) 
initial  radius  of  mixed  region 
mean  temperature 
time  to  start  of  collapse 
mean  velocity  vector 
velocity  defect 

velocity  defect  at  the  centerline 
body  velocity 

streamwise,  horizontal  cross-stream  and  vertical 
cross-stream  mean  velocity  components 


fluctuating  velocity  vector 


u,v,v 


streamwise,  horizontal  cross-stream  and  vertical 
cross-stream  fluctuating  velocity  component 
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u,v,w  root  mean  square  fluctuating  velocity  components 

W wake  width 

x*  position  vector 

x,y,z  streamwise,  horizontal  cross-stream  and  vertical 

cross-stream  coordinates,  respectively 

Z wake  half-width 


Zl/2 


distance  between  centerline  and  point  where  k = k(0)/4 


Greek  Symbols 
a 
6 
e 
0 


v ,v 

V W 


V .V 

Tv*  Tw 


P 


temperature  gradient  of  undisturbed  fluid 
volumetric  expansion  coefficient 
turbulent  energy  dissipation  rate  per  unit  volume 
temperature  fluctuation 

eddy  diffusivity  coefficients  for  momentum  in  the 
horizontal  and  vertical  direction,  respectively 

eddy  diffusivity  coefficients  for  heat  in  the 
horizontal  and  vertical  directions,  respectively 

density 

density  at  the  centerline 
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INTRODUCTION 

A turbulent  wake  in  a stably  stratified  fluid  exhibits  some  peculiar 
features  which  are  not  encountered  in  wakes  in  nonstratified  media. 

Whereas  axisymmetric  wakes  in  homogeneous  fluids  increase  indefinitely 
in  width  under  the  action  of  the  diffusive  turbulent  forces,  wakes  in 
stably  stratified  environments  are  known  to  reach  a maximum  vertical 
extent  and  subsequently  collapse  to  a minimum  height.  At  the  same  time, 
under  the  combined  action  of  turbulent  diffusion  and  gravitational  forces, 
they  increase  in  the  horizontal  direction  faster  than  they  would  in  a 
neutral  medium.  Furthermore,  since  the  decrease  in  the  vertical  extent 
is  associated  with  a reduction  in  potential  energy,  gravitational  internal 
waves  are  generated,  which  propagate  through  the  surrounding  fluid  (See 
Fig.  1.) 

The  first  experimental  investigation  of  this  phenor,  ...  m ■?  a s co"  ced 
by  Schooley  & Stewart  [1],  who  used  a flow  visualizaf* . ?::•  nniqu?  to 

observe  the  wake  behind  a self-propelled  body.  The  ScJu  fix ..  baj  since 
been  documented  in  detail  by  Stockhausen,  Clark  & Kennedy  [2],  Lin  & Pao 
[3] , [4]  and  Stromm  [5],  whereas  Pao  & Lin  [6],  [7]  focused  on  the  wake 
behind  a towed  body  in  a stably  stratified  medium.  Van  der  Watering,  Tulin 
& Wu  [8],  Sumdaram  et  al.  [9]  and  Meritt  [10]  investigated  the  collapse 
of  mixed  turbulent  regions,  generated  by  different  means  in  quiescent, 
stably  stratified  media.  These  phenomena  are  closely  related  to  the  decay 
of  momentumless  wakes  in  such  media,  since  turbulence  in  a momentumless 
wake  is  known  to  behave  in  a way  similar  to  turbulence  behind  a point 
source.  Wu  [11]  has  observed  the  collapse  of  a tube  of  homogeneous,  non- 
turbulent  fluid  in  a stably  stratified  environment  and  the  resulting  wave 


I 


i 

j 

* 

i 

1 


i 


pattern. 
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Theoretical  investigations  by  Wessel  [12]  and  Young  & Hirt  [13]  focus 
on  the  collapse  of  a nonturbulent  mixed  region;  this  flow,  however, 
differs  from  the  collapse  of  the  turbulent  mixed  region:  whereas  the 
latter  reaches  a final  finite  width,  the  former  collapses  indefinitely, 
until  the  fluid  in  the  mixed  region  reaches  the  equilibrium  level.  The 
only  theoretical  treatment  of  the  mixed  region  is  by  Vasiliev  et  al.  [14]. 
The  turbulence  model  used  by  these  authors,  however,  does  not  predict  any 
collapse,  but  only  a slowing  down  of  the  rate  of  growth  of  the  turbulent 
region.  Furthermore,  they  neglect  the  nonlinear  terms  of  the  Navier 
Stokes  equations.  It  is  not  clear  that  this  approximation  is  justified 
because  these  are  the  very  terms  that  cause  the  collapse.  Finally,  they 
make  no  comparison  with  experimental  data  and  give  no  predictions  for  the 
streamwise  velocity  component  in  wake  flows. 

In  this  work  it  is  proposed  to  use  a simple  turbulent  flow  model, 
coupled  with  a computational  procedure,  which  will  enable  one  to  describe 
and  document  the  phenomenon  of  collapse  of  turbulent  wakes  in  stratified 
media. 

THE  FUNDAMENTAL  EQUATIONS 

Consider  the  wake  behind  a body  moving  with  uniform  velocity  in  a 
stably  stratified  medium  with  uniform  temperature  (or  density)  gradient. 
(The  application  of  the  model  to  salinity-stratified  media  is  straight- 
forward). Assuming  the  Boussinesq  approximation  to  be  valid  (i.e.,  that 
the  specific  weight  of  the  fluid  is  variable  but  its  density  can  be 
considered  constant)  one  obtains  the  following  form  for  the  Reynolds 
stress  and  continuity  equations 
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In  addition,  the  following  simplifying  assumptions  will  be  made 


(a)  The  flow  Is  elliptic  In  z and  y and  parabolic  In  x 
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By  differentiating  Eq.  (6)  with  respect  to  y and  Eq.  (7)  with  respect 
to  z and  adding,  a Poisson  equation  for  the  pressure  can  be  obtained: 


£+''2>+  J£+»2>--68iT<T--)  • 


In  addition,  the  following  equation  governs  the  transport  of  temperature: 


„ 3T  . 3T  . „ 3T  3 -r-  3 -5- 

°o  3^  + vW  + w • 37  ve  - a?"6 


In  order  to  solve  the  above  system  of  equations,  closure  approximations 
are  needed  for  the  Reynolds  stresses  uv  and  uw,  as  well  as  for  the  turbulent 
heat  fluxes  v0  and  w0.  These  will  be  discussed  in  the  following  section. 


THE  TURBULENT  FLOW  MODEL 

Recently,  the  Reynolds  stress  model  has  become  increasingly  popular 
in  computational  investigations  of  turbulent  flow.  For  flows  in  wakes 
behind  both  towed  and  self-propelled  bodies,  however,  the  simpler  k-e 
model  gives  results  which  are  in  better  agreement  with  experiment.  (See 
Hassid  [15]).  Needless  to  say,  of  course,  the  k-e  model  is  also  preferable 
because  of  its  lower  requirements  in  computer  storage  and  time. 

The  closure  used  here  is  essentially  based  on  the  model  of  HcGuirk 
and  Rodi  [16].  Transport  equations  are  used  for  both  the  turbulent  energy 
k and  the  dissipation  rate  e and  gradient  diffusion  forms  for  the  Reynolds 
stresses.  However,  the  expressions  for  the  latter  quantities  are  derived 
from  invariant  second  order  models  by  neglecting  the  advection  and 
turbulent  transport  terms. 

The  expressions  for  u^Uj  and  u^0  are 


McGuirk  and  Rodi  used  their  model  to  describe  jet  flows  in  stratified 
media.  Wake  flows  differ  from  jet  flows  in  that  the  latter  are  parabolic 
in  y.  Therefore,  in  our  case,  some  further  approximations  are  needed. 

As  already  stated,  vw  will  be  neglected  in  comparison  to  the  normal  Reynolds 
stress  components.  Thus,  the  equations  for  the  turbulent  heat  fluxes  are 
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These  equations  will  be  further  simplified  by  assuming  that  3T/3y 
is  much  smaller  than  3T/3z.  Although  this  assumption  is  not  strictly  true, 
the  error  caused  will  be  small,  both  for  weak  stratification,  because 
the  dependence  on  3T/3y  is  small  in  that  case,  and  for  strong  stratifi- 
cation, because  then  3T/3z  is  indeed  much  larger  than  3T/3y.  Thus  Eq.  (14) 
becomes 


cT  e 3z 


Substituting  Eq.  (18)  into  Eq.  (17)  one  finds 


(1“C2T)  k2  „ 3T 

i+^r~78^ 

~2  ~~1 

Now,  from  Eq.  (11),  one  obtains  the  following  expressions  for  v and  w“ 


2 2 (l-ci~c2^ 

v ’ ' 3 “ Oj  k > 


w 3 c, 


2 = - 7= — — k - 2 6 g - w6  . 


Substituting  into  Equations  (16)  and  (19)  one  obtains 


-E  „ 3T  k 3T 

- * ■ V 37  ■ c T 37  • 


vQ  - v 

Tw  3z 


id  ii 

c8  e 3z 

, . . k2  a 3T 
1 + A^  6 8 3l 


. * 


f.l£  -f  iifch 
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where 


and 


2 (l-crc2) 
'0  3 Cj^  cT 


2 ^~C2T  ^~c2 

A - — [ — + -] 

C1T  CT  C1 


(24) 


(25) 


Turning  to  the  equations  for  the  Reynolds  stress  one  has 


uv  “ - 


1-c  2 

2 k v 3U 

Ci  e 3y 


(26) 


uw  ■ - 


(^•— CJ«)  . 2 «„  

[!UL|2 . 2 8 g u9] 


(27) 


Again,  using  Equations  (20)  and  (26)  one  obtains 


where 


— „ 3U  k2  3U 

~UV“  vv  37“  ^7-97  * 


2 (l-c2)  (l-c^c.,) 

®M  " " 3 cl  cl 


(28) 


(29) 


Neglecting,  as  already  mentioned,  3T/3y  in  the  equation  for  u0  and 
substituting  in  Eq.  (27),  one  ends  with 


where 


and 


— 3U  k_  (14A*F)  3U 

UW  w 3z  = Si  e (1+A"F) (1+AF)  3z 


F - 6 E — — 
P 8 2 3z 

e 

1-c 


A’  - A - (2 


2 + 2T.  1 

C1  C1T  C1T 


A" 


1-c, 


C1T  C1 


(30) 

(31) 

(32a) 

(32b) 
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Here  the  model  differs  from  the  one  of  McGulrk  and  Rodl  because  these 
authors  claim  that  gravitational  production  does  not  appear  In  the  equation 
for  e.  It  was  found  that  such  an  assumption  leads  to  unrealistically  low 
values  for  k.  It  should  be  mentioned  that  there  is  wide  a-  greement 
among  different  research  groups  on  the  question  of  the  influence  of 
buoyant  forces  on  the  e equation. 

The  constants  used  "re 


I 


|| 

I 

I 1 


Of  all  the  approximations  used  the  neglect  of  3T/9y,  which  was 
necessary  in  order  to  obtain  the  gradient  forms  for  the  Reynolds  stress 
equations,  is  the  one  likely  to  cause  the  greatest  error.  Notice  however, 
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that  for  the  non-buoyant  case,  the  expressions  derived  above  are  exact. 

METHOD  OF  SOLUTION 

The  system  of  equations  to  be  solved  consists  of  the  x and  y momentum 
equations 


3Ud  3Ud  9Ud  3 ,Ps  2.  3 .2  2. 

Uo37-  + V3r  + W3r-  + ^(F+v)  + 9^(u  "V) 


a a 9U^ 

f-  [v  3-^]  + [v  ■—} 

3y  v 3y  3z  w 3z 


(35) 


„ 3V  . „ 3V  . „ 3V  . 3 ,Ps  . 2.  _ 

U -s—  + V^—  + W-r—  + -s-  (—  + v ) *=  0 , 

o 3x  3y  3z  oy  p 


(36) 


the  transport  equation  for  temperature 


U il+v  — + W — = — (V  51)  + 1_  (V  W) 

o 3x  V 3y  W 3z  3y  ^ Tv  3y'  3z  1 Tw  3z' 


(37) 


the  Poisson  equation  for  pressure 


P a J\2  P ry  r\2  9 9 ^ 

” (tt  + v'  > + <—  + wZ)  + *-=■  (w  -v  ) = B g 4-  (T-az)  , (38) 

3yZ  P 3zZ  P 3zZ  dZ 


the  continuity  equation  (Eq.  8),  the  transport  equations  for  the  turbulent 

energy  k (Eq.  33),  and  the  dissipation  rate  t (Eq.  34).  In  addition,  the 

2 2 2 2 

following  expressions  for  u -v  and  w -v  can  be  derived  from  Eq.  (11) 


and 


2 2 (1-c2)  . — 3U  — 3UX 

u — (uv  — + uw  ^) 


2 2 — 
w — v = 


(39) 


B g w8 


(40) 
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On  AB, 


3W  3 2.  9Ud 

v ■ 37 ' 37  (r  + v > ' IT  ' 0 ■ 


On  AD, 


3V  3 ,Ps  2.  SUd  • - 

37-w-3i(r  + ”)  = 3r'°  • 


On  BC, 


V - -p  - P - U . = 0 . 

tfz  s d 


On  CD, 


This  system  was  converted  into  a set  of  difference  equations.  Forward 
differencing  was  used  in  the  x direction,  centered  differencing  for  the 
diffusive  terms  and  upwind  differencing  for  the  lateral  transport  terms. 

The  mesh  system  used  (Fig.  2)  consists  of  an  inner  mesh  AEFG  in 
which  the  turbulent  processes  occur  and  an  outer  mesh,  in  which  there 
is  no  turbulence.  The  outer  mesh  was  three  times  as  large  as  the  inner 
mesh.  The  reason  is  that  both  experimental  evidence  and  computer 
calculations  suggest  that  the  waves  caused  by  the  collapse  of  a non turbulent 
mixed  region  advance  twice  the  original  width  of  that  region  after  two 
Brunt-Vaisala  periods.  Since  there  are  very  few  experimental  data  for 
times  larger  than  two  Brunt-Vaisala  periods  and  in  addition,  the  wave 
phenomena  are  much  weaker  when  the  mixed  region  is  turbulent,  this  mesh 
was  estimated  to  be  adequate. 

The  mesh  used  was  allowed  to  expand,  so  that  the  original  growth 
due  to  turbulent  diffusion  can  be  simulated. 

The  following  boundary  conditions  were  used: 


For  k and  c,  the  boundary  conditions  were 
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On  EF  and  FB,  k ■ e ■ 0 . (45) 

On  *E,  ||  - § - 0 . (46) 

On  AG,  fj-lf-0  • (47) 

Admittedly,  some  of  these  conditions  are  artificial  and  may  result  in 
errors  as  the  wave  fronts  approach  the  mesh  boundaries.  It  is  estimated, 
however,  that  these  errors  will  not  be  excessive  for  times  less  than  two 
Brunt-Vaisala  periods.  As  a check,  the  differential  equations  were  solved 
using  different  boundary  conditions  and  this  did  not  significantly  affect 
the  results  in  the  turbulent  region. 

RESULTS  AND  COMPARISON  WITH  EXPERIMENT 

The  main  scaling  parameter  used  by  most  investigators  (Refs.  3,  6,  8, 

9 and  10)  for  wake  flows  in  stratified  media  is  the  Froude  number, 
defined  as 


Fr 


(48) 


where  N is  the  Brunt-Vaisala  frequency 


N 


1_ 

2tt 


/*!£)!/ 2 
(P  3yJ 


(49) 


However,  this  scaling  factor  is  not  sufficient;  what  determines  the 
extent  that  a momentumless  wake  will  grow  before  collapsing  is  rather  the 
Richardson  number,  which  expresses  the  ratio  between  the  potential  energy 
and  the  turbulent  kinetic  energy  at  the  origin 
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The  decay  of  momentumless  wakes,  even  in  the  self-similar  state, 
depends  on  such  factors  as  the  body  shape  (whether  the  body  is  blunt  or 
slender),  the  propelling  mechanism,  etc.  For  non-buoyant  flows  behind 
slender  bodies,  the  wake  width  is  given  by 


H p *x. • 22 

D * 8 V 


(See  Lin  & Fao  [3]).  In  Fig.  3a  the  influence  of  stratification  on  the 


wake  behind  such  a body  is  shown.  It  is  seen  that  the  wake  height  and 

width  become  functions  of  Nx/U  . Collapse  starts  when  Nx/U  = 1/3.  The 

o o 

wake  height  decreases  by  10%  and  then  increases  slightly.  The  wake  width 
increases  throughout  the  flow. 

In  the  same  figure,  the  correlations  of  Lin  & Pao  [3]  are  shown.  The 
agreement  between  the  predicted  and  the  measured  wake  width  is  very  good; 
however,  our  predictions  overestimate  the  height  of  the  wake.  It  must  be 
noted,  however,  that  Lin  & Pao  concede  that  there  is  considerable  difficulty 
in  estimating  the  wake  width  because  of  the  irregular  fluctuating  nature 
of  the  outer  wake  boundary. 

In  Fig.  3b  the  width  and  the  height  of  the  wake  behind  a towed  body 
in  stratified  flow  are  shown.  Here  the  computed  values  of  2z^y2  are  compared 
with  the  values  of  H and  W measured  by  Pao  & Lin  [5]  who  claim  that  for 
drag  wakes  H “ ^z\f 2'  our  conPutations>  21  = .)  It  is  seen  that 

in  wakes  the  first  maximum  is  reached  at  a higher  value  of  Nt  than  in 
momentumless  wakes.  Though  the  computations  fail  to  reproduce  the  10% 
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reduction  In  width  recorded  In  the  experiments,  the  agreement  can  be 
considered  satisfactory,  in  view  of  the  influence  of  stratification  on 
the  initial  wake  widths,  which  cannot  be  accounted  for  by  the  calculation 
procedure  and  in  view  of  the  experimental  errors  of  flow  visualization. 

In  Figs.  4a  and  4b,  the  predicted  wake  width  is  compared  with  the 
results  of  Lin  & Pao  [4]  for  momentumless  wakes  and  Pao  & Lin  [7]  for 
drag  wakes.  The  system  of  maxima  and  minima  observed  in  these  experiments 
is  reproduced;  there  is  however,  a certain  difference  in  phase  probably 
due  to  the  use  of  V=0  and  W=0  as  initial  conditions. 

In  Fig.  5a  and  5b  it  is  seen  that  stratification  results  in  slower 
decay  of  the  velocity  defect,  both  for  drag  wakes  the  momentumless  wakes. 
Agreement  with  the  experimental  data  is  very  good. 

In  Figs.  6a  and  6b  the  decay  of  turbulent  energy  with  and  without 
stratification  is  shown.  It  is  seen  that  in  spite  of  the  additional 
destruction  of  turbulent  energy  due  to  the  stable  stratification  the 
difference  in  the  turbulence  levels  is  very  small,  especially  in  the 
momenturaless  wake.  Agreement  between  theory  and  experiment  is  again 
very  good  in  spite  of  the  fact  that  the  predicted  ^ is  compared  with  the 
measured  u. 

Finally  in  Table  I the  predicted  values  of  the  maximum  and  the 
minimum  heights  as  well  as  the  collapse  times  of  mixed  regions  and 
momentumless  wakes  are  compared  with  their  values  as  measured  by  several 
experimental  workers. 

In  this  table,  tcQ^  is  the  time  (for  mixed  regions;  distance  divided 

by  mean  body  velocity  in  momentumless  wakes)  at  which  the  wake  reaches 

the  first  maximum  height,  Z and  Z . are  the  maximum  and  minimum  wake 

max  min 

heights,  respectively,  whereas  the  Richardson  number  R is  defined  as 
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ro  ^8  ag 
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(52) 


Table  1 


Collapse  Time,  Maximum  Height  and  Minimum  Height  of 
Turbulent  Mixed  Regions  in  Stably  Stratified  Media 


Reference 

'co/6 

Z 

max 

Zmin 

Meas. 

Pred. 

Meas. 

Pred. 

Meas. 

Pred. 

Schooley  & Stewart 

Fr  «* 

60 

2.12 

1.8 

2.16 

2.00 

1.36 

1.8 

Stockhausen,  Clark 

& Kennedy 

Fr  = 

100 

1.63 

1.2 

3.00 

2.08 

2.10 

1.8 

Van  der  Watering, 

1 /? 

Tulin  & Wu 

R1/2 

= 1.34 

1.33 

.92 

1.51 

1.26 

1.11 

.96 

1.56 

1.48 

.98 

1.25 

1.23 

1.04 

.91 

1.56 

1.21 

.98 

1.31 

1.26 

1.06 

.93 

2.19 

1.35 

.90 

1.11 

1.17 

.90 

.80 

2.19 

1.35 

.91 

1.19 

1.17 

. .90 

.82 

3.09 

1.40 

.78 

1.06 

1.09 

— 

.68 

4.23 

1.23 

.51 

.95 

1.03 

.77 

.57 

4.27 

1.05 

.69 

1.00 

1.05 

.83 

.62 

1.56 

1.12 

.93 

1.24 

1.20 

.94 

.86 

i 

1.57 

1.12 

.94 

1.25 

1.21 

.89 

.87 

2.28 

.89 

.74 

1.20 

1.13 

1.00 

.72 

3.21 

1.25 

.93 

1.13 

1.13 

.94 

.63 

2.06 

1.05 

.65 

1.20 

1.06 

1.03 

.63 

• 

1.89 

.73 

.70 

. 1.04 

1.07 

.88 

.65 

3.69 

1.07 

.50 

1.04 

1.09 

— 

.58 

5.21 

1.28 

.58 

.88 

1.04 

.77 

.59 

5.29 

1.18 

.45 

.95 

1.02 

.88 

.55 

Sundaram,  Straton 

3.20 

3.75 

.72 

1.54 

1.12 

1.07 

.71 

& Rehm 

— 

2.80 

2.40 

.77 

1.31 

1.13 

1.00 

.74 

where  Vq  is  the  initial  rate  of  spreading  and  rQ  is  the  initial  radius  of 
the  mixed  region.  For  the  experiments  of  Schooley  & Stewart  and  Stockhausen, 


Clark  & Kennedy,  rQ  is  the  body  diameter. 
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There  are  of  course  several  sources  of  error  in  this  table.  Since 
no  initial  turbulent  energy  data  were  available,  this  quantity  was  deduced 
from  the  initial  rate  of  spread  and  the  equation 

— XT  = *26  (53) 

which  is  applicable  to  self-similar  flow  of  momentumless  wakes  in 
nonstratified  media.  The  initial  turbulent  energy  and  dissipation  profiles 
were  also  assumed  self-similar.  Given,  however,  the  very  high  Richardson 
numbers,  it  is  very  unlikely  that  the  turbulent  energy  distribution  may 
reach  self  similarity.  Besides,  the  predicted  value  of  tco^  is  the  value 
of  t when  the  wake  height  actually  starts  decreasing;  typically,  however, 
it  has  decreased  by  less  than  5 Z at  t=2  tcQ^  and  it  is  only  then  that  it 
can  be  observed  in  a flow  visualization  experiment.  Finally,  for  the 
experiments  of  Schooley  & Stewart  and  Stockhausen,  Clark  & Kennedy,  the 
production  of  turbulence  by  mean  velocity  gradients  could  not  be  taken  into 
account. 

In  spite  of  all  these  considerations,  it  is  thought  that  in  view 
of  the  crudeness  of  the  observation  methods,  the  numerous  assumptions  needed 
for  the  calculation  and  the  disagreement  between  the  data  themselves,  even 
the  order  of  magnitude  agreement  presented  in  Table  1 is  satisfactory. 

DISCUSSION  AND  CONCLUSIONS 

The  development  of  momentumless  wakes,  drag  wakes  and  turbulent  regions 
in  stably  stratified  media  was  investigated  using  a simple  turbulent 
energy  - dissipation  model.  This  is  the  simplest  model  that  can  be  used 
for  such  flows;  evidence  from  calculations  in  nonstratified  media  suggests 
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that  for  wake  flows,  this  model  Is  also  more  accurate  than  the  more 
sophisticated  Reynolds  stress  models. 

The  effect  of  stratification  is  expressed  through  a damping  factor 
in  the  eddy  diffusivity  coefficients  for  vertical  exchange  of  momentum 
and  heat.  The  eddy  diffusivities  in  the  horizontal  direction  remain 
unaffected. 

The  stabilizing  temperature  gradient  acts  in  two  ways:  (a)  by  reducing 
the  turbulent  fluctuation  level  in  the  vertical  direction,  thus  impeding 
the  ability  of  turbulence  to  diffuse  in  this  direction;  (b)  by  generating 
a mean  flow  pattern  in  which  the  vertical  velocity  is  directed  towards 
the  horizontal  symmetry  plane  and  the  horizontal  velocity  is  directed  away 
from  the  vertical  symmetry  plane.  This  flow  pattern  is  due  to  the  difference 
in  weight  between  the  fluid  in  the  turbulent  region  and  the  surroundings. 
Thus,  expansion  of  the  mixed  regions  is  impeded  in  the  vertical  direction 
and  enhanced  in  the  horizontal  direction.  The  wake  eventually  reaches  a 
maximum  height  in  the  vertical  direction  and  shrinks.  A system  of  maxima 
and  minima  in  the  height  follows,  which  are  separated  by  roughly  one 
Brunt-Vaisala  period. 

The  stratification  is  shown  to  slow  the  decay  in  the  velocity  defect, 
presumably  as  a result  of  the  decrease  in  the  rate  of  vertical  exchange 
of  momentum.  It  has  very  little  effect  on  the  total  turbulent  energy, 
because  the  destruction  of  turbulent  energy  by  gravitational  forces  is 
balanced  by  the  increase  in  the  production  (since  the  rate  of  decay  of  the 
mean  velocity  is  smaller)  and  the  decrease  in  the  diffusion.  These 
features  of  the  flow,  which  are  well  documented  in  the  experiments  of  Lin 
& Pao,  are  well  reproduced  by  the  model. 
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The  model  has,  however,  several  limitations.  It  is  not  known  whether 
it  can  account  for  the  laminarization,  that  Lin  & Pao  claim  occurs  after 
four  Brunt-Vaisala  periods.  The  influence  of  gravitational  negative 
production  of  turbulent  energy  on  the  dissipation  equation  is  not  certain 
and  the  particular  assumption  used  in  this  model  is  only  a tentative  one, 
to  be  modified  when  a more  realistic  model  is  available.  The  calculation 
procedure  used  is  not  reliable  for  times  larger  than  two  Brunt-Vaisala 
periods,  because  at  such  times  the  artificially  imposed  boundary  conditions 
may  result  in  reflection  of  the  waves  generated  by  the  collapse.  In  order 
to  overcome  this  drawback,  one  may  introduce  a region  of  high  viscosity 
on  the  outer  part  of  the  flow,  so  that  outcoming  waves  are  damped  before 
they  can  be  reflected.  If  this  is  done  one  may  have  a more  reliable 
view  of  the  wave  pattern  caused  by  the  collapse  but  the  computer  require- 
ments in  storage  and  time  would  be  significantly  increased.  Finally 
one  should  mention  the  artificial  viscosity  errors  caused  by  the  upwind 
differencing  technique  used  to  approximate  the  lateral  transport  terms 
which  may  partially  account  for  the  fact  that  the  predicted  collapse  rate 
is  lower  (i.e.,  the  final  thickness  is  higher)  than  in  many  (but  not 
all)  experiments. 
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(b)  Drag  Wake,  Data  of  Pao  & Lin  [7] 

Figure  4 - Height  of  Wakes  in  Stratified  and  Nonstratified  Media 
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